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ABSTRACT 

One remaining open question regarding the physical properties of Lya emitters (LAEs) is their dust 
content and its evolution with redshift. The variety of results is large and with those reported by 
now is difficult to establish clear relations between dust, other fundamental parameters of galaxies 
(star-formation rate, metallicity or age) and redshift. In this Letter, we report Herschel PACS-100/mi, 
PACS-160/im and Spitzer MIPS-24/im detections of a sample of spectroscopically GALEX-selected 
LAEs at z~0.3 and ~1.0. Five out of ten and one out of two LAEs are detected in, at least, one 
PACS band at z~0.3 and ~1.0, respectively. These measurements have a great importance given that 
they allow us to quantify, for the first time, the dust content of LAEs from direct FIR observations. 
MIPS-24/im detections allow us to determine IR properties of the PACS-undetected LAEs. We obtain 
that mid-IR/FIR detected star-forming (SF) LAEs at z~0.3 have dust content within 0.75;$ ^ 1200 ^ 

<2.0, with a median value of A x ~1.1. This range broadens out to 0.75 < A ° <2.5 when 
~ ' 1200A to ~ 1200A ~ 

considering those LAEs at z^l.O. Only one SF LAE is undetected both in MIPS-24/im and PACS, 

with A ° 5:0.75. These results seem to be larger than those reported for high-redshift LAEs and, 
1200A ^ 

therefore, although an evolutionary trend is not clearly seen, it could point out that low-redshift LAEs 
are dustier than high-redshift ones. However, the diverse methods used could introduce a systematic 
offset in the results. 

Subject headings: galaxies: evolution — galaxies: stellar content — infrared: galaxies — ultraviolet: 
galaxies 



1. INTRODUCTION 

Lya emitters (LAEs) have been proved to be of great 
importance in the study of the formation and evolution of 
galaxies through cosmic time. At 2<z<7, a large number 
of candidates have been found, mainly via the narrow 
or intermediate band technique in the optical /NIR 
wavelength range (lOuchi e t al. 2008; Cowie & Hu 1998; 
Gronwall et al. 2007; Guaita et all [20 10: Shi oya et all 



20091: iM uravam a et al.l 120071: Nilsson et all 12009: 
Fuiita et al.l 12003: B ongiovanni et al.l l2010l ). At z<2, 
the LAEs catalogued up to now have been found 
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by spectroscopic s urveys using the Galaxy Evolution 
Explorer (GALEX, iMartin et all 12005( 1. The selection 
technique is based on looking for a Lya line in the 
UV spectra of objects with a measured U V continuum 
(jDeharveng et al.ll2008l : ICowie et al.ll2010D . 

Despite of the large number of LAE candidates found 
at different redshift s, their physical properties and their 
evolution with cosmic time are not completely under- 
stood yet. Moreover, the relation between LAEs and 
other kind of high redshift galaxies, such as Lyman break 
galaxies (LBGs) or distant red galaxies (DRGs), remains 
unclear. Some studies have suggested that the stellar 
populations of LAEs exhibit a double nature, some be- 
i ng old and massive and other youn g and less massiv e 
(jGawiser et al.l [20061 : lLai et al.l 12003 lOno et al.ll20Toh . 
However, recent studies have proved that stellar popu- 
lations of LAEs vary continuously and smoothly over a 
wide range, confo rming a big zo o of galaxies (Oteo et al. 
2011, submitted; iNilsson et all 120 111 ). 

One remaining open question is whether LAEs are 
primitive and dust-free objects, evolved and dusty, or 
present a wide range of dust attenuation. The Lya line 
has a resonant nature and, therefore, the probability of 
being absorbed by the interstellar medium (ISM) even 
before leaving out the host galaxy is quite high. Due 
to that, it was initially thought that LAEs found via a 
narrow band technique should be young and dust free 
objects, allowing the L ya line to appear in emission. 
However, Neufeld (|1991l ) proposed a scenario where the 
Lya equivalent width could be enhanced by the pres- 
ence of dust when it resides in cold and neutral clumps 



2 



Oteo et al. 



located within an ionize d ISM. Actually, s ome studies 
(jFinkelstein et alJl2QQ9al l3 iBlanc et alJl201Q[ ) have obser- 
vationally confirmed this hypothesis, opening the pos- 
sibility of finding the Lya emission in non primordial 
galaxies at different redshifts. The dust content in LAEs 
has been studied in a wide range of redshifts and dif- 
ferent results have been found, from dust-free objects to 
low, moderate, and even high dust content. In Table [TJ 
some results on the study of dust attenuation in LAEs 
at different redshifts are summarized. 

Most previous works have studied the dust attenuation 
in LAEs by fitting the UV and optical rest-frame spectral 
energy distrib ution (SEP) with st e llar p opulation tem- 
plates, mainly iBruzual fc Charlotl (|2QQ3|) (BCQ3 ) ones, 
and assuming the dust attenuation law of lCalzettil (p.997). 
However, both dust and old stars can redden the rest- 
frame integrated UV continuum, varying the expected 
colors from the dust-free scenario. UV and optical pho- 
tons which are absorbed by dust are reemitted in the 
FIR region. Therefore, the most accurate way to study 
the dust properties of high red shift galaxies is o b serving 
the dust emissio n in th e FIR. iFinkelstein et al.l (j2009b) 
and iDayal et"aTI (j2Q10h studied the expected detection 
of dust emission from LAEs at z>4 and z~6, respec- 
tively, finding that a large percentage of them at those 
redshifts should be detected in only a few hours with 
ALMA. IFinkelstein et al.l (j2QQ9bh also found that LBGs 
are 60% more likely to be detected than LAEs, meaning 
that they are dustier objects and showing a possible dif- 
ference between the two populations. However, some re- 
cent studies did not find FIR counterparts for LBGs z~3, 
chara cterizing their FIR SEP only with stacked ph otom- 
etry dMagdis et al.l l2Q10b iRigoDoulou et "all l2010h. For 
LAEs. iBongiovanni et al.l (j2010f ) reported the detection 
in the FIR of 3 L AEs at z~2 .3 sele cted from the AL- 
HAMBRA survey (jMoles et al.ll2QQ8h . being their opti- 
cal to mid-IR SEP w ell fitted by AGN-like templates of 
iPolletta et all (|2QQ7D . Oteo et al. (2011a, submitted), 
also found FIR counterparts for 5 spectroscopically con- 
firmed LAEs at z~2.5. In this case, according to their op- 
tical spectra, only one is an AGN, and the other four are 
star-forming (SF) LAEs. All of these FIR-detected LAEs 
resulted to be (U)LIRGs candidates and, therefore, red 
and dusty objects, despite having the L ya line in emis- 
sion a nd with large equivalent widths. 1 Chapman et al.l 
(2005f) also found evidences of the Lya emission in a sam- 
ple of sub-mm galaxies detected in 850/im, being most 
of them (U)LIRGs candidates too. 

Until now, few FIR counterparts for LAEs at z<2 have 
been reported, and for those LAEs with FIR counter- 
parts, an extensive analysis of their FIR properties has 
not been done. In this letter, we report the FIR de- 
tections of a sample of LAEs at z~0.3 and z~1.0 by 
using data coming from Photodetector Array Camera 
and Spectrometer (PACS. IPoglitsch et al. 2010) on board 
Herschel Space Observatory. These FIR measurements, 
combined with UV information, enable us to determine 
their dust content in an accurate way without the uncer- 
tainties coming from the properties of the dust itself, as 
its emission in the FIR, missing in the optical methods. 

This Letter is organized as follows: In Sections [2] and 
[3] we present the LAE sample and the FIR data used, re- 
spectively. In Section 3] we report the FIR counterparts 
of our LAEs and in Section [5] we study their dust attenu- 



ation, comparing them with previous results. In Section 
[6j the main conclusions of our work are described. 

Throughout this work we assume a flat universe with 
(O m ,OA,^o) = (0.3,0.7 ,0.7), and all mag nitudes are 
listed in the AB system (jOke fc Gunnlll983l ). 

2. LAE SAMPLE 

In our analysi s , we use an LAE sample selected 
from iCowie et all (j2010[ h where the spectra, taken with 
GALEX NUV and FUV grisms, were selec ted from the 
Multimission Archive at STScI (MAST). Cowie et al. 
(|2010f ) used an automatic search procedure to look for 
emission lines in the spectra of objects with a measured 
UV continuum. For each source, a polynomial curve was 
used for constraining the continuum, and a gaussian pro- 
file fit to look for emission lines. Given that the spec- 
tra become very noisy at the edges of the spectral range 
in each channel, only those objects which have the Lya 
emission in the wavelength ranges 1452.5-1750 A and 
2006-2735 A for the FUV and NUV grisms, respectively, 
were selected. In the redshift space, this means that the 
sources are located within z=0. 195-0.44 (z~0.3 sample) 
and z=0.65-1.25 (z-1.0 sample) in the FUV and NUV 
spectra, respectively. Both samples are nearly complete 
to an NUV limiting magnitude of 21.5 and have Lya lu- 
minosity limits of 5-10 41 and 1-10 43 erg cm _1 s _1 at z~0.3 
and ~1.0, respectively. LAEs were classified as AGN or 
star-forming (SF) galaxies according to the width and 
shape of the Lya line and the presence or absence of 
AGN ionization lines in the UV spectra. Figure [T] shows 
the spatial distribution of GALEX-selected LAEs in both 
COSMOS and GOODS-South fields. 

3. FIR DATA 

GOODS-South and COSMOS fields have been ob- 
served with PACS on board Herschel in the framework 
of PACS Evolutionary Probe project (PEP, PI D. Lutz). 
PEP is the Herschel Guaranteed Time Key-Project to 
obtain the best profit for studying FIR galaxies evolu- 
tion from Herschel instrumentation (Lutz et al. 2011, in 
prep). For both fields, PACS-100/im and PACS-160/im 
observations are available, with 3a limiting fluxes of (1.1, 
5.0) in PACS-100/im and (2.0, 11.0) in PACS-160/im in 
(GOODS-South, COSMOS). PACS-100/im and PACS- 
160/im bands cover the FIR wavelength ranges 80-130/im 
and 130-2 10/im, respectively. Two kind of catalogues 
have been bu ilt: blind and with MI PS-24/i m position pri- 
ors by using iLe Floc'h et al.l ()20Q9h and iMagnelli et all 
(|201lh catalogs in COSMOS and GOODS-South, respec- 
tively. In order to have as much FIR information as 
possible and more reliability in the identification of the 
sources, we use the catalogs extracted with MIPS-24/im 
priors. Data reduction, catalogs construction and simu- 
lations, aimed at deriving completeness, fraction of spu- 
rious sources and photometric reliability, are described 
in Lutz et al. (2011, in prep.). 

Figure [T] also shows the PACS- Herschel sources ex- 
tracted with 24/im positions priors. It can be seen that 
PACS observations in COSMOS field cover the whole re- 
gion where GALEX-selected LAEs are located, whereas 
the surveyed area by PACS in GOODS-South only cov- 
ers a little portion of GALEX observations. In this 
sense, 28 and 4 GALEX-selected LAEs are within the 
PACS surveyed area in COSMOS and GOODS-South, 
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respectively. Table [2] summarizes the number of GALEX- 
selected LAEs within the PACS area in both fields. 

4. LAES WITH FIR COUNTERPARTS 

In order to find out how many LAEs have FIR coun- 
terparts we performed a match between GALEX-LAEs 
coordinates and PACS catalogs in both fields. We look 
for a possible FIR detection within a radius of 2" around 
the position of each LAE, since it is the astrometric un- 
certainty in the position of the sources. With this proce- 
dure, we find that a noticeable fraction of LAEs, mostly 
above 30%, have PACS-FIR detections in at least one 
considered PACS band, considering limiting fluxes shown 
in Sect. [U In Table [2j the number of PACS detection 
are also shown, separated by field, redshift and nature 
of the object (AGN or SF). The percentages of FIR de- 
tections are larger in GOODS-South than in COSMOS 
because of the depth of the observations. If the PACS 
observation had been in COSMOS as deep as in GOODS- 
South, a large number of FIR counterparts would have 
been expected. 

A source confusion analysis is needed to ensure that 
each FIR detection corresponds to the matched GALEX- 
selected LAE. We take 30" optical cut-outs around the 
position of each GALEX-selected LAE and overplot the 
MIPS-24/im contours in order to verify if the detections 
in each band is related or not to the same object. We 
use the MIPS-24/im contours given that the PACS cata- 
logues are based in 24/mi positions priors, where the PSF 
of the observations is ~30% narrower than redder bands 
and, therefore, the accuracy in the positions is better. It 
can be seen in Figure [2] that each FIR detection reported 
corresponds to the matched GALEX-LAE. 

5. DUST ATTENUATION IN STAR-FORMING LAES 

The ratio between the total IR and UV luminosities, 
Lir and L[/y, respectively, is directly related to the dust 
attenuation in galaxies, and calibrations have been a l- 
ready built (jMeurer et al.lll999l : iBuat et al.ll2005l . I2007D . 
Since all the PACS-detected SF GALEX-LAEs also have 
UV measurements, we can infer their dust attenuation 
from direct observations of their FIR SED. We use the 
NUV band given that, on one hand, it is not contami- 
nated by the Lya emission for those GALEX-LAEs at 
z~0.3 and, on the other hand, it is not affected by the 
Lyman break for those GALEX-LAEs at z^l.O. \jjjv 
can be obtained directly from the observe d magnitudes 
in the NUV band, as in IBuat et al.l (|2005f ). L/# can be 
calculated by using the calibrations given in Oteo et al. 
(2011a, submitted) between PACS bands and total in- 
frared luminosities: 

log L IR = 0.99 log L 100fim + (0.44 ± 0.25) (1) 

\ogL IR = 0.96 log L 160/im + (0.77 ± 0.21) (2) 

where all the luminosities are in solar units and Lioo/xm 
and Li60/im are defined as vL v . PACS-detected LAEs 
have measurements in 100/im and/or 160/im. For those 
detected in both bands, we obtain the Lm by using the 
calibration in the PACS- 160/im band, since it is the clos- 
est one to the dust emission peak at our redshift ranges 
and, therefore, produces a more accurate result. 



We can also obtain the total IR luminosities, i.e. dust 
content, for those PACS-undetected but MIPS-24/mi de- 
tected SF GALEX-LAEs by using the Chary & Elbaz 
(2001) calibrations between the total IR and 15/im and 
12/im luminosities for local galaxies. These wavelengths 
are within the rest-frame ranges sampled through the 
MIPS-24/im band for LAEs at z-0.3 and -1.0, respec- 
tively. The applicability of the se calibrations to high- 
redshift galaxies is discussed by Elbaz ^et al.l ((2010) who 
found that, below z—1.5, mid-IR extrapolations to the 
total IR luminosities are correct for SF galaxies below 
the (U)LIRG regime. The MIPS-24/im in formation was 
taken from the from S-COSMOS survey ([Sanders et al.l 
2007). Only one object is not detected in that band. 

In Figure [3] we plot the ratio Lir/Luv ( or ? equiva- 
lently, dust attenuation) versus Lir of our mid-IR/FIR 
detec ted GALEX-LAEs, b esides the corresponding data 
of ([Gil de Paz et al.1l2007l hereafter G07) for a collection 
of nearby, SF galaxies. Our LAEs are located in the 
low-dust regime of the G07 data, indicating that objects 
exhibiting the Lya line in emission at low and moderate 
redshift should not be as dusty as the whole population 
of galaxies at the same redshift. The dust attenuation in 
our mid-IR/FIR-detected LAEs at z~0.3 ranges within 
0.75<A x 5:1-8 mag, with a median value of 1.1 mag. 

~ 1200A^ to 

This range broadens up to A s 5:2.5 when consider- 
to ^ 1200A ~ 

ing the z-1.0 LAEs as well. For the GALEX-LAE un- 
detected both in MIPS-24/im and PACS we assume an 
upper limit in its dust attenuation equal to the minimum 
value obtained for the mid-IR/FIR detected GALEX- 
LAEs: A * < 0.75. 
1 i^noA 1 , 1 

IFinkelstein et al.l (|2011[ ) studied dust content in LAEs 

at z~0.3 by using SED fitting and Balmer decrement. 

They derived an attenuation in the stellar continuum 

from Balmer decrement ranging from A 2 ~0 to ~1.8, 

to to 1200A ' 

with a median value of 0.78 mag. These values are 
in agreement with those found via SED fitting, be- 
ing most of them w ithin the combined la uncertain- 
ties. iCowie et al.l (|2010l ) obtained dust attenuation 
from the NUV continuum slope and Balmer decrement 
of their LAEs at z~0.3, obtaining a median value of 
A i2oo j^~2.25 and ~1.8, respectively, which are notice- 
ably high^r__th^Ji__th^ in Finkels tein et al.l 
(|2011f ). iNilsson fc Mollerl (|2009f ) reported significantly 

higher values, most of them being A 2 >2.0, and as 
to > t> 12 ooA^ ' 

large as 7.5 mag. The results derived here are s lightl y 

higher tha n those deri v ed in IFinkelstein et al.l ((2011), 

similar to Cowie et al.l (120101) v alues, and lower than 

those in Nils son fc M0llerl (j2009[ ). These discrepancies 

may be caused by the different methodologies employed, 

and even when they are same, the results do not agree 

quite well, likely due to the low number of objects in each 

sample. 

The dust attenuation values reported up to z~2 are 
similar to those obtained in this work. However, at z>3, 
while there are some measurements of high dust attenu- 
ation in high redshift LAEs, the vast majority of the re- 
ported extinction measurements are substantially lower 
than A 2 50.75 and are, in many cases, consistent 

1200A^ ' J ' 

with essentially no extinction. This could indicate that 
high redshift LAEs are less dusty than those at lower 
redshift s. However, note that, overall at the highest red- 
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shifts, the number of the studied LAEs at each redshift 
is low and, therefore, the results are not statistically sig- 
nificant. Furthermore, it is worth noting that SED fit- 
ting estimations (sometimes done with stacking analy- 
sis), spectroscopic results and IR measurements could 
produce systematic differences. This means that a pos- 
sible evolution of dust attenuation with redshift should 
be taken with caution. 

6. CONCLUSIONS 

In this letter, we have reported the FIR detections of 
a sample of GALEX-selected LAEs at z~0.3 and z^l.O. 
Our main conclusions are as follows: 

• We find that a large fraction of LAEs at z~0.3 and 
z^l.O are detected in the FIR. In the worst case, 
the fraction of FIR detected LAEs is about 15%, 
although in most cases the fractions are higher than 
30%. Therefore, up to z^l.O, at least 30% of the 
LAEs are expected to be detected in the FIR under 
the depth of the catalogs used. The percentages of 
MIPS-24//m detections are even higher, only one 
object being undetected in that band. 

• Based on direct measurements of the mid-IR/FIR 

SED, we conclude that the studied SF-LAEs at 

z~0.3 present a wide range of dust attenuation, 

0.75<A i <1.8, with median value A j of 
~ 1200A^ ' 1200A 

1.1. This range broadens up to 0.75<A 2 ^2.2 

to 1 ~ 1200A^ 

when considering LAEs at z^l.O. We only find one 

object, undetected both in MIPS-24/im and PACS, 

whose dust content has an upper limit of A 2 < 

11 1200A ~ 

0.75. 

• The dust attenuation values derived here are com- 
parable to those reported in previous works up to 
z~2.0, pointing out that there is no significant evo- 
lution in dust content in that epoch. At z>3, al- 
though some works showed the existence of dusty 
LAEs, most of the extinction measurements are 



lower than those at z<2.0, which could indicate 
an evolution with redshift. However, the differ- 
ent methods used to derive the results and the low 
number of LAEs at certain redshifts make the re- 
sults statistically unsignificant. 
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Fig. 1. — Spatial distribution maps of the GALEX-LAEs (both AGNs and SF) and PACS sources. Blue diamonds and blue dots are 
PACS-undetected and PACS-detected GALEX-LAEs, respectively. Red dots are PACS detections with 24/xm priors. Left: COSMOS field. 
Right: GOODS-South fields. It can be seen that all the GALEX-LAEs in COSMOS are within the covered area of PACS. In GOODS- 
South, only four GALEX-LAEs are within the PACS imprint. Note in the left plot that one of the PACS-undetected sources is so close to 
a PACS-detected one that the point is not clearly seen. 




Fig. 2. — 30" x 30" ACS-optical cut-outs and MIPS-24/i flux isocontours, linearly scaled to the maxima. GALEX-LAEs are the objects 
located in the center of each image. Black curves are flux iso-contours in MIPS-24/im band, linearly scaled to the peak flux. 
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Fig. 3. — Ratio between total infrared and NUV luminosities against total i nfrare d luminosity for star-forming GALEX-LAEs. Right 
vertical axis (dust atten uation) was built by using the calibration of Buat et al. (2005). Yellow dots are data of nearby galaxies taken from 
IGil de Paz et al. ( 2007). Dark shaded zone is the average region (zbl.5cr) where these objects lay. Light shaded region corresponds to the 
dust content range of LAEs at z~0.3 reported by Finkelstein et al. (2009a). Red and blue dots are SF LAEs at z~0.3 and z~1.0, respectively. 
Filled symbols are PACS-dectected LAEs and open symbols are PACS-undetected LAEs. We do not plot the MIPS-24 j um-undetected and 
PACS-undetected SF-LAE at z~0.3 (see text) given that there is no way to obtain its total IR luminosity with the calibrations used. 
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TABLE 1 

Summary of some previous results in dust attenuation in LAEs. Reference, redshift, dust attenuation range and the 
method followed to obtain it are shown: sed fittings with bc03 templates to the individual or stacked observed 
photometry (TF-BC03 and TF-BC03 & Stacked, respectively), optical spectroscopic analysis (Spec), UV continuum slope 

AND MID-IR MEASUREMENTS. 



Reference 



redshift Dust attenuation [mag] 



Method 



Finkelstein et al. (2009a) 
Finkelstein et al. (2011) 
Nilsson fc Moller (2009) 
Cowie et al. (2010) 
Cowie et al. (2010) 
Guaita et al. (2011) 
Nilsson et al. (2011) 
Nilsson fe Moller (2009) 
Gawiser et al. (2007) 
Lai et al. (2008) 
Ono et al. (2010) 
Ono et al. (2010) 
Nilsson et al. (2007) 
Ono et al. (2010) 
Ono et al. (2010) 
Finkelstein et al. (2008) 
Finkelstein et al. (2009c) 
Pirzkal et al. (2007) 
Lai et al. (2007) 
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Spec 

IRAC-8/im measurements 
UV continuum slope 
Spec 

TF-BC03 

TF-BC03 
MIPS-24/im measurements 
TF-BC03 & Stacked 
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TF-BC03 & Stacked 
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TABLE 2 

Number of GALEX-selected LAEs within the surveyed area of PACS in COSMOS and GOODS-South fields. We separate 

THE OBJECTS IN REDSHIFT, NATURE (ACCORDING TO THE UV SPECTRA, SEE TEXT) AND THE FIELD WHERE THEY ARE LOCATED. NUMBERS 

BETWEEN PARENTHESES ARE THOSE DETECTED IN THE FIR WITH PACS. 



Field 


Type 


z-0.3 


z-1.0 


COSMOS 


AGN 


3(1) 


14 (2) 


COSMOS 


SF 


9(4) 


2(1) 


GOODS-South 


AGN 


0(0) 


3(2) 


GOODS-South 


SF 


1 (1) 


0(0) 



